Shear instabilities in granular mixtures 
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Dynamical instabilities in fluid mechanics are responsible of a variety of important common phe- 
nomena, such as waves on the sea surface or Taylor vorteces in Couette flow. In granular media 
dynamical instabilities has just begun to be discovered. Here we show by means of molecular dynam- 
ics simulation the existence of a new dynamical instability of a granular mixture under oscillating 
horizontal shear, which leads to the formation of a striped pattern where the components are seg- 
regated. We investigate the properties of such a Kelvin-Helmholtz like instability and show how it 
is connected to pattern formation in granular flow and segregation. 
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A simple example of instabilities in fluid mechanics is 
the Kelvin-Helmholtz instability where a flat interface 
between two fluids flowing one past the other at different 
velocities is unstable, explaining for instance why waves 
on the see surface form . In athermal systems such as 
granular mixtures, where an hydrodynamic like theory is 
not yet well established, the discovery of instabilities is a 
recent achievement [j. Q S By studying the evolution 
of two bands of different grains placed on an horizontally 
oscillating tray (see Fig. 1(a)), we discover a new dynam- 
ical instability. This instability appears as growing wavy 
interface, which, at long times, leads to a pattern of alter- 
nating segregated stripes of grains, perpendicular to the 
driving direction. We suggest that this type of surface 
shear instability is the common mechanism for appar- 
ently different phenomena such wave formation as found 
in experiments on granular flow Iff, and segregation pro- 
cesses 7] as observed in 0, 0, lid ll l| , initially explained 
in terms of a thermodynamic driven phase separation. In 
our simulations we determine the region of existence of 
the pattern formation - segregation process as a function 
of the area fractions of large and small grains, and we 
discuss the dependence of the pattern as a function of 
the amplitude and frequency of oscillation. 

MD simulation model - We perform Molecular Dy- 
namics (MD) simulations (also known as Discrete Ele- 
ment Methods) of a binary system of disks which lay on 
a tray (i.e., in two dimensions) with periodic boundary 
conditions in the x-direction and hard walls in the other. 
Grains interact with the tray by a viscous force propor- 
tional to their relative velocity via a viscosity parameter 
/jl different for the two species; when overlapping, grains 
interact via 2D Hertzian contact forces 12]. Simulations 
details and parameters are chosen to model experimen- 
tal conditions similar to those investigated by Mullin and 
coworkers HH El El 
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2D this reduces to the linear spring-dashpot model, f n = 
k n SijUij — 7n^redV n ij, where k n and j n are the elastic 
and viscoelastic constants, and m re d = m^mj/ (m^ + rrij) 
is the reduced mass. As in EH we model the interaction 
with the tray via a viscous force f t = — /i(v— v tray ), where 
Vtray(^) = ^ttAu sin(z4)x is the velocity of the tray and 
v the velocity of the disk, plus a white noise force £(£) 
with (£(£)£(£')) = 2r <K* - t'), where T = 0.2 g 2 cm 2 S - 3 . 
The interaction between particles and walls is elastic. We 
solve the equations of motion by Verlet algorithm with 
an time-step of 6/is. For the grain-grain interaction, we 
use the value k n = 2 10 5 g cm 2 s -2 and 7 n chosen, for 
each kind of grains, such that the restitution coefficient 
is given: e = 0.8 14]. The two components of our mix- 
ture (named b and s) have viscous coefficient /i& = 0.28 g 
s _1 and n s = 0.34 g s _1 . Apart from a simple rescaling 
of masses and lengths, these values are those of Ref. E5 
(and given in private communications), and are taken 
from direct measurements on the experimental system 
in |g. The size of the tray is L x = 320 cm, L y — 16 
cm. The qualitative picture we discuss does not change 
if these values are changed. 

Results - We consider by now a mixture of heavy 
grains, of mass M5 = lg and area fraction fa = 0.37, 
and light grains, with M s = 0.03 g and (j) s = 0.41, all 
with the same diameter D = 1 cm, prepared in a hori- 
zontally fully segregated configuration (Fig. 1(a)). The 
mixture is vibrated on an horizontal tray along the x- 
direction with amplitude A and frequency v. The char- 
acteristic time scale, 75 = Mb /fib and r s = M s //i s , is 
different for the two species and they are forced to os- 
cillate with different amplitudes and phases, subject to 
an effective surface shear periodically varying in time. 
One may expect configurations in which the two differ- 
ent components form two stripes parallel to the driving 
direction (see Fig. 1(a)) to be stable, the two stripes 
oscillating independently with different amplitude and 
phases. Instead the initially flat interface between the 
two components is observed to evolve via the formation 
of a surface modulation which has a growing amplitude. 
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FIG. 1: (color on-line) Evolution of a binary mixture of heavy (red-dark gray) and light (blue-ligth gray) grains on an 
horizontal tray oscillating along the x-direction with amplitude A = 1.2 cm and frequency v — 12 Hz (only a quarter of the 
system length is shown in each figure). Panels a The system is initially prepared in a horizontally segregated configuration 
with a fiat interface between the two species. Under shaking it develops a surface instability with a sine like modulation of 
growing amplitude leading at long times to a state where segregated stripes perpendicular to the driving direction appear. 
Upper small panels show the corresponding evolution of the surface instability model Eq.s JIJ in presence of an oscillating 
shear Av(t) = 2-kAv sin(27rz/£). Panels b The same system is now initially prepared in a mixed state and then shaken. The 
formation of local fluctuations of the density of the two species (shown at t = 1.6 s) leads to microsurfaces generating the same 
instability mechanism as in panels a. 



Finally, this leads to a state where segregated stripes ap- 
pear, perpendicular to the driving direction (at variance 
with known results in colloidal fluids This insta- 

bility appears of the same kind of that observed in two 
fluid systems and in liquid sand systems subject to hori- 
zontal oscillations under gravity, which is responsible for 
the ripples observed on the shoreline In these cases, 
however, because of gravity there is an energy cost as- 
sociated to the growth of the interface, which therefore 
gets stabilized. 

When the system starts from a disordered initial state, 
as shown in Fig. 1(b), the instability develops too. The 
formation of local fluctuations of the density of the two 
species create microsurfaces which are forced to evolve 
according to the instability shown in Fig. 1(a). For in- 
stance Fig. 1(b), at t= 1.6 s, shows the existence of 
clusters of particles of the same species elongated per- 
pendicularly to the driving direction. As time goes on 
these clusters merge giving rise to the striped pattern of 
Fig. 1(b) (t= 80 s). It is apparent that the originally 
disordered system is ordered, i.e., segregated, by the dy- 
namic instability. Fig. 1(a) (t= 320 s) and Fig. 1(b) 
(t= 80 s) show that the wavelength of the final steady 
state depends on the initial conditions. 

The segregation process shown in Fig . 1(a) reproduces 
that observed experimentally in || ^| . In the present 
set-up grains radii are equal, no "depletion" forces can 
be at work, and no "thermodynamic" coarsening phe- 
nomena can be considered responsible for this segrega- 
tion process, as initially conjectured [ll| by making an 
analogy with colloidal systems. Such a remark appears 
to be consistent with a scenario proposed by use of mod- 
ified Navier- Stokes equations Interestingly, segrega- 
tion from a disordered initial state to a final striped one 



has also been observed in a two fluid system oscillated in 
absence of gravity [l^ . 

The above dynamical instabilities should manifest in 
the mixture also under different shearing protocols. To 
check this, we study now a different configuration where 
gravity induces granular flow and thus a differential shear 
on the mixture components appears. The tray no longer 
oscillates, but it is inclined of an angle with the respect 
to the horizontal by rotation around the y-axis.The two 
species of grains attain different limit velocities (of the 
order of ~ gr^ sin and v s ~ gr s sin 9, with g gravi- 
tational acceleration) and thus experience a differential 
shear. As shown in Fig. 2, under these conditions the 
initially flat interface between the mixture components 
changes too, via Kelvin-Helmholtz-like instabilities. An 
unstable pattern forms and curl-like structures appear. 
This phenomenon, induced by surface dynamical insta- 
bilities, is very similar to those recently discovered in 
experiments on granular flows |6|. 

The results of Fig. 1 and Fig. 2 point out the essen- 
tial role of shear induced surface instabilities in pattern 
formation and the connections with segregation transi- 
tions. Even though this is not an usual fluid, the com- 
plex nature of this phenomenon can be understood via 
schematic hydrodynamic-like considerations. On the two 
sides of the perturbed interface flux lines are narrowed 
and widened as the pressure is respectively lowered and 
increased (Bernoulli law). The pressure gradient rein- 
forces the perturbation of the interface giving rise to a 
positive feedback. Within such a perspective, the essen- 
tial features of pattern formation observed in the two 
experiments of Fig.s 1 and 2 are captured by a model 
originally proposed for the Kelvin-Helmholtz instability 
of granular flows |6|- The model describes the evolution 
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FIG. 2: (color on-line) Evolving interface in the same granular mixture of Fig. la flowing down a tray rotated of an angle 
= 12° with respect to the horizontal around the y-axis. The red (dark gray) grains are faster that the blue (light gray) 
ones. Upper small panels The upper panels show the evolution, at corresponding times, of a Kelvin-Helmholtz-like surface 
dynamical instability modeled by Eq.s (QJ in presence of a constant shear Av(t) — Avo, reproducing a pattern qualitatively 
similar to the one found in the main panels. 



of the velocity components, vx and vy, of the system 
interface initially located at Y(X) = 0: 

v x = Av(t) tanh(F) v Y = csin(k x X) (1) 

where Av(t) is the relative bulk velocity between grains 
far on the two side of the interface at time t, as kx and 
c are constant parameters (we consider the same val- 
ues used in 0|, i.e., kx = 5, c = 0.1, but our results 
are robust to changes). The presence of a velocity gra- 
dient, Av(t), independently from its microscopic origin 
(different grains friction, different driving of the com- 
ponents, etc.), enhances shear and results in patterns 
formation. In the case of gravity induced deformations, 
Av(t) = Avq = Vb — v s can be considered to be time 
independent since in the stationary regime the two limit 
velocities are given. In such a case, the model essen- 
tially coincides with the one introduced in 6] and well 
describes the qualitative features of the MD simulations 
of Fig. 2 (see upper panels). In the case of the vibrated 
horizontal tray, the relative bulk velocity is oscillating in 
time and we fix Av(t) = 2ttAv sin(27rz4). This gives rise 
(see upper panels in Fig. 1(a)) to a growing interface with 
features very close to those observed in the simulations. 

Pattern formation - We show now under which condi- 
tions pattern formations occur and explain how the prop- 
erties of the final state of the mixture, such as the emerg- 
ing characteristic length scale of the stripes, A, depend 
on both the dynamics control parameters and the rel- 
ative concentration of the two components. This gives 
rise to a complex mixing/segregation diagram which we 
discuss in details here in the case of the horizontally 
vibrated experiment. In order to make direct compar- 
isons with experimental results on segregation observed 
by Mullin and coworkers || EJ [HI , we consider now a 
binary mixture of large mono-disperse disks, of diameter 
Db = 1cm, covering an area fraction </>5, and small poly- 
disperse disks |2l|, with average diameter D s = 0.7cm 
(with 17% polydispersity) and area fraction <fr 3 . Masses 
and frictional parameters are as given before. 

The system starts from a mixed initial configuration 
as in Fig. 2(b). The nature of the state reached by the 



system under shaking at late stage is crucially depen- 
dent on the area fraction of the two species: segregation 
in stripes is only found for high enough concentrations. 
This behaviour, in the (0&,0 S ) plane, is summarized in 
the diagram of Fig. 3(a) showing the system "fluid" and 
"crystal" regions along with their segregation properties, 
for v — 12 Hz and A — 1.2 cm. Large grains are con- 
sidered to be in a "fluid" configuration when their radial 
density distribution function, g(r), shows a first peak at 
r = Db and a second one at r = 2D5, and to be in a 
"crystal" configuration when a new peak at r = V3Db 
appears . The system is in a "glassy" state ji^ when 
on the longest of our observation time scales, the system 
is still far from stationarity. 

Fig. 3(a) shows that grains at small concentrations 
are mixed and in a fluid state. Segregation via stripes 
formation appears at higher concentrations. At even 
higher concentrations, large grains form stripes with a 
crystalline order, as smaller grains are always fluid for 
their polydispersity. Finally, at very high area fractions, 
the system is blocked in its starting disordered config- 
uration ("glassy" region). For instance, by increasing 
(j) s at a fixed value of (\>b (say <pb 0.174), we observe 
first a transition from a mixed fluid state to a segregated 
striped fluid and then a transition where the the monodis- 
perse phase crystallize. The experiments of [ll|, where 
(j)b c± 0.174, show the very same transitions found here at 
locations differing by a 10%. 

In the case <pb = 0.30 and (j) s = 0.28, where stripes 
form, we describe their dependence on the dynamics con- 
trol parameters in Fig.s 3(b) and 3(c), showing that the 
length scale, A, increases as a function of the shaking fre- 
quency, and of the amplitude A. These results are to 
be compared, for instance, with those found in liquid- 
sand mixtures under oscillating flow, where ripples form 
with a wavelength depending on the amplitude of oscilla- 
tion, but not on its frequency [2^. The dependence on v 
can be schematically understood by comparison with the 
characteristic time scales T5 = Mb /fib and r s = M s /fi s of 
the two species (here r 6 _1 = 0.28 Hz and r s _1 = 11.3 Hz): 
in the limit v ^> t^ -1 ,?-" 1 grains are not able to follow 
the tray motion and no sensitivity to v is expected, as 
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FIG. 3: Panel a Ordering properties of the late stage config- 
urations of the mixture as a function of the area fractions of 
the two components. The shaded area covers the region where 
segregation via stripes formation occurs, circles: large grains 
are in a fluid state, squares: large grains form a crystal, 
stars: the system appears blocked in a "glassy" disordered 
configuration (see text). When stripes form their character- 
istic length scale A is a function of the frequency, i/, and of 
the amplitude, A, of the driving oscillations. This is shown, 
in the case (j)b — 0.30 and (j) s — 0.28, in panels b and c. 

well as when v <C r" 1 since the grains move with the 
tray. Analogously, the dependence on A is expected to 
be substantial when A is at least of the order of the mean 
grains separation length, I = (Afa/irDl + 4</> s /7i\D;?) _1//2 , 
since under this condition grains strongly interact. 

In conclusion our simulations of a binary mixture un- 
der horizontal vibrations have revealed the existence of a 
new dynamical instability and shed light on the process 
of size segregation under oscillatory shear and its con- 
nections to pattern formation in granular flows. Within 
such a unifying framework, we derived the "phase di- 
agram" of the mixing/segregation states of the mixture 
and its corresponding transitions. Finally, such hydrody- 
namic like processes appear to be related to those known 
in thermal fluids and fluid-grains systems, ranging from 
Kelvin-Helmholtz instabilities to ripples formation in liq- 
uid sand mixtures, even though interesting difference are 
found. 
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